SNe observations in a meatball universe with a local void 



Kimmo KainulainerQ and Valerio MarrgfO 

Department of Physics, University of Jyvaskyla, PL 35 (YFL), FIN-40014 Jyvaskyla, Finland and 
Helsinki Institute of Physics, University of Helsinki, PL 64, FIN-00014 Helsinki, Finland 

We study the impact of cosmic inhomogeneities on the interpretation of observations. We build 
an inhomogeneous universe model without dark energy that can confront supernova data and yet 
is reasonably well compatible with the Copernican Principle. Our model combines a relatively 
small local void, that gives apparent acceleration at low redshifts, with a meatball model that gives 
sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different 
effects of voids on the data, allow the model to mimic the concordance model. 
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Introduction. The "safe" consequence of the success 
of the concordance model is that the flat isotropic and 
homogeneous ACDM model is a good phenomenological 
fit to the real inhomogeneous universe. We will discuss 
the possibility that the late time evidence for dark energy 
can be explained by the coincidental late time formation 
of nonlinear large scale inhomogeneities. Indeed, until 
the nature of the dark energy is completely understood, 
it is useful to look for alternative models that fit the data. 

We will assume that the spacetime of the inhomoge- 
neous universe is accurately described by small perturba- 
tions around the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) solution whose energy content and spatial cur- 
vature are denned as Hubble- volume spatial averages over 
the inhomogeneous universe. Following Ref. pQ the lat- 
ter will be called the Global Background Solution (GBS). 
The cosmological model obtained through observations, 
on the other hand, will be called the Phenomenological 
Background Solution (PBS). The idea is to associate the 
concordance model with the PBS, while the actual GBS 
is the Einstein-de Sitter (EdS) model. 

Setup. For consistency with the CMB spectrum and 
the age of the universe, the current background expan- 
sion rate of the EdS model is taken to be = 100 
km s _1 Mpc -1 with = 0.5. The local expansion rate 
Hq will be higher, about h ~ 0.58 in the current model, 
which is within 2-cr of the HST key project [2] value of 
0.72 ±0.08 [3 j. Our setup is made of two elements. First, 
we will model the overall universe by a meatball model [4 
consisting of two families of randomly placed halos. The 
family A describes very large, low density contrast struc- 
tures and the family B models large clusters of galaxies. 
The parameters specifying the model are (for each fam- 
ily) the average comoving distance between meatballs d c , 
the proper radius of the meatball R p and the mass of the 
meatball M. The numerical values we used are given in 
Table H| The distance d c between meatballs is connected 
to their comoving number density by n c = j^f^yd~ 3 . 
We assume the meatballs to have virialized at z = 1.6 
and therefore their proper radii R p are constants. For 
the family B we used the singular isothermal spheres 



TABLE I: Parameters of the meatball model (ao = 1). 



Quantity 


Family A 


Family B 


d c 


100 h' 1 Mpc 


10 h' 1 Mpc 


R p 


10 h' 1 Mpc 


580 h' 1 kpc 


M 


6.1 • KF/i^M© 


6.1 • lO^h^Me 


density profile 


Gaussian 


SIS 



(SIS) density profile and defined the radius to have a 
density contrast of 200 at virialization. For the family 
A the gaussian profile was used and the radius was cho- 
sen to qualitatively describe the largest structures seen 
in simulations of structure formation. The distribution 
of family A meatballs leaves underdense regions of order 
100/i _1 Mpc which are filled by the family B meatballs. 
Because the mass is equally subdivided between the two 
families, these underdense regions have on average a den- 
sity contrast of 5 « —0.5. The second element of our 
setup consists of placing the observer in one of the large 
underdense regions and modelling the metric of this par- 
ticular void around the observer more accurately with 
a Lemaitre-Tolman-Bondi (LTB) bubble matched to the 
EdS background metric. See Fig. [I] for a sketch. 

Lensing. In a late-time universe dominated by voids 
the homogeneity is recovered only on scales larger than 
the largest inhomogeneity scale, which in our model is 
Lhom ~ 100/i _1 Mpc. However, type la supernovae (SNe) 
probe angular scales LsNe <^ ^hom an d it is therefore not 
clear if the physics inferred from SNe observations can be 
directly associated with the smoothed-out GBS model, 
especially with small data samples. Epigrammatically, 
the commutativity between averaging and measuring is 
not guaranteed [lj. Note that this idea, which defines 
the weak-backreaction, is different from the possible non- 
commutativity between averaging and dynamics which is 
the kernel of strong-backreaction studies [5]. 

Cumulative gravitational lensing is one possible source 
for the non-commutativity between averaging and mea- 
suring and to study this effect we have focused on a 
meatball model in this letter. This model incorporates 



FIG. 1: The larger disks represent the meatballs of the family 

A, while the smaller ones represent the meatballs of the family 

B. The concentric circles mark the overdense shells in the 
LTB bubble and are roughly 100/i _1 Mpc apart. 



quantitatively the crucial feature that photons can travel 
through voids and miss the localised over densities. This 
feature, instead, is not present in swiss-cheese models 
where the bubble boundaries are designed to have com- 
pensating overdensities. Such models have indeed been 
shown to have on average little lensing effects [3 [7] . In 
the weak-lensing approximation the lens convergence k 
is given by 



k(z) = I ^ } drG(r : r s (z))5(r : t(r)), 
Jo 

where <5(r, t) is the density contrast and 
3H^ r(r s - r) 1 



G(r,r s ) 



2c 2 r s a(t(r)) ' 



(1) 



(2) 



The functions a(t), t(r) and r(z) correspond to the EdS 
background model, r s is the comoving position of the 
source at redshift z and the integral is evaluated along 
the unperturbed light path. Neglecting the second-order 
contribution of the shear, the shift in the distance mod- 
ulus caused by lensing is 



Am(z) = 51og 10 (l - k(z)) . 



(3) 



Am is proportional to the total matter column density 
along the photon path. For a lower-than-GBS column 
density the light is demagnified (empty beam S — — 1, 
for example), while in the opposite case it is magnified. 

In Ref. [8] we derived a fast and easy way to obtain 
the convergence probability distribution function (PDF) 
for meatball models. First define a normalized profile 
(p(x) = p(x)/p, where p is the density profile of a meat- 
ball and p is the EdS density. Then the following quan- 
tity characterizes the contribution from one meatball, hit 
with an impact parameter 6, to the convergence: 



r(M) 



R(t) 



2xdx 



\[x* 



b 2 



(4) 



FIG. 2: Top: the 307 SNe of the Union Compilation of Ref. [9] 
binned with Az = 0.1. Bottom: shift in the distance modulus 
Am for No taken from the top panel and for the meatball 
model of this paper. 



Then divide the comoving distance r s to the source and 
the radius R of the meatball into Ns and Nr bins of 
widths Ari and A6 m , such that Ab m <C R and R <C 
Ari <C r s . The convergence due to a meatball placed 
within the bin (i, m) is just 



«iim = G(ri,r a )r(b m ,ti) . 



(5) 



The PDF is then generated by the following functional 

N s N R ,^ v 

K({hm}) =yzyz n ^™> ( -^p- - AN im ) , (6) 

where ki m is a Poisson random variable of parameter 
Nq ANim. No(z) is the number of observations at red- 
shift z and AN im = n c AV irn is the expected number of 
meatballs in the volume AV irn = 27rb rn Ab rn Ari. For a 
Poisson variable of parameter A, the mean is A and so 
the expected convergence in Eq. ([6| is zero, consistent 
with photon conservation in weak lensing. The PDF is 
generated by creating a large sample of configurations 
{ki m } drawn from the Poisson distribution, and comput- 
ing the appropriate convergences through Eq. This 
is a task that a laptop fullfils in less than a second as 
opposed to expensive ray tracing techniques. The weak 
lensing approximation of Eq. ([6| describes the mode of 
the PDF very well always and, in case of non-point-like 
meatballs, also the tails of the PDF with < 5% of error. 

Eq. ^ displays explicitly the effect of the size of the 
data sample: even if n might have a skewed PDF and a 
nonzero mode for Nq = 1, the gaussianity is recovered for 
large No and the set of observations becomes unbiased. 
To define No(z) we will use the Union Compilation of 
Ref. [9 , which consists of 307 SNe spread between < 
z < 1.6. We have binned the SNe with a bin width of 
Az = 0.1 and the result is plotted in the top panel of 
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Spherically symmetric LTB metrics can be written 



as 



FIG. 3: AH r /H (top panel) and density contrast (bottom 
panel) with respect to normalized proper distance (R/a(t)) 
at times corresponding to the redshifts indicated in the plot. 



Fig. [2] SNe observations within the same bin are then 
be treated as repetition of the same observation. We 
have used a modification of the turboGL package [10] to 
evaluate the magnification bias for our model and for the 
Union Compilation. Results are shown in the bottom 
panel of Fig. |2) The magnification bias is nonneglible for 
z > 0.5 and large for z > 1. We stress that the latter 
bias is statistical and not coming from selection effects. 

Hubble bubble. The second element of our model con- 
sists of using an LTB bubble matched to the EdS back- 
ground of the meatball model to describe the local metric 
around the observer. The new ingredient is that an ob- 
server inside a void expanding faster than the background 
sees an apparent acceleration (see, e.g., Refs. [TTJ [12]). 
This effect is easy to understand: our cosmological ob- 
servables are confined to the light cone and hence tempo- 
ral changes can be associated with spatial changes along 
photon geodesies. For example "faster expansion now 
than before" is simply replaced by "faster expansion here 
than there". This is why a local hubble bubble model 
can mimic the effect of a cosmological constant at re- 
cent times (z < 1). The price to pay in a simple hubble 
bubble model is that the inhomogeneity has to extend 
to the point in space/time where the effect of dark en- 
ergy disappears. This requires [12] an enormous local 
void of radius 1.5 — 2 /i _1 Gpc. Moreover, to avoid a too 
large dipole in the CMB, our position in the void would 
have to be very special leading to a gross violation of the 
Copernican principle. The model of this letter features 
inhomogeneities on much smaller scales. Here the effects 
of the dark energy at large redshift (z > 0.5) are mim- 
icked by the statistical lensing bias, and the local void 
only needs to model the cosmological constant at small 
redshifts where the lensing effects are neglible. 



ds 2 = _ c 2^2 + R' 2 {rt) + r2 
1 — k(r)r z 



dr 2 +R 2 (r,t)dn 2 . (7) 



This reduces to the usual FLRW metric when R(r, t)/r = 
a(r,t) — > a(t) and k(r) — > ±1 or 0. In the dust case the 
Einstein equations for this metric give 



a 2 {r,t) 
a 2 {r,t) 



8nG a/ , c 2 k(r) 



(8) 



where p is the average density up to the shell r. We 
specify our LTB model at the recombination time with a 
density contrast of order ~ 10 -3 in accordance with the 
CMB. In particular, the curvature function is defined by 
2 overall parameters that specify size and depth of the 
local void and by 3 parameters that specify the detailed 
profile inside the void which is shown by the evolution 
of the density contrast in Fig. [3] Also shown is the ra- 
tio AH r /H where AH r = H r - H, H is the EdS ex- 
pansion rate and H r = R! / R' . At late times the initial 
perturbations generate two fast expanding (AH r > 0) 
voids of radius 50 — 100 /i _1 Mpc, surrounded by two thin 
collapsing (AH r < 0) shells. The inner void of radius 
^ 100/i _1 Mpc, where the observer is located, has a local 
expansion rate of h ~ 0.58 which, as anticipated, is larger 
than the EdS value. Moreover, the density contrast of the 
latter is close to S ~ —0.5 consistently with our meatball 
model for the overall universe. In order not to have a 
too large dipole in the CMB the observer has to be at a 
distance < 15 /i _1 Mpc from the center of the inner void. 
Tighter constrains will come from future probes sensitive 
to off-center anisotropy (see, for example, [13 ). 

Results. Let us now describe the main features of our 
results shown in Fig. [4] First, since our observer is inside 
a negative-curvature dominated region, the slope of the 
inhomogeneous model starts almost flat at z = 0, as it 
is evident in the zoom of the bottom panel. However, at 
finite redshifts the spatial variation of the expansion rate 
mimics the apparent acceleration and a positive Am pro- 
portional to the difference between local and background 
expansion rates Ho — Hqo appears. Beyond the dimen- 
sions of the local hubble bubble, at z > 0.1, the slope 
of the curve follows the EdS model until the lensing bias 
effect becomes important at z > 0.5 after which it starts 
to pull the curve up again as shown in the top panel of 
fig. [4] This smooth rise of the predicted Am at large z 
and the behaviour around z ~ 0.1 — 0.2, also clearly vis- 
ible in Figure [4j are the main characteristic differences 
between our model and the ACDM model. Both these 
issues would be easily resolved by a JDEM-like survey 
with a data set of 2000 SNe. 

We also computed the x 2 m order to estimate the 
goodness-of-fit of our model. It should be stressed that 
we did not explore the parameter space to minimize the 
X 2 but chose the parameters essentially by hand. Thus 
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FIG. 4: Top panel: shown is the distance modulus with re- 
spect to the empty universe for the inhomogeneous model of 
this paper, for the ACDM model and for the EdS model to- 
gether with the full Union Compilation of Ref. 9 . Bottom 
panel, zoom for low redshifts. 



the comparision between our model and the concordance 
model is qualitative only. Numerical results are shown in 
Table [Til Obviously the simple EdS model has a poor fit, 
but introducing inhomogeneities improves the x 2 signifi- 
cantly. The local Hubble bubble has a larger impact on 
X 2 , but lensing also has a large positive effect. 

TABLE II: X 2 for Union Compilation of 307 SNe. 

Model X 2 

ACDM 312 

ACDM + meatballs 323 

EdS 608 

EdS + H. bubble 440 

EdS + H. bubble + meatballs 396 



Finally, including inhomogeneities changes also the 
predictions of the ACDM model. The second line in Ta- 
ble HIl shows results of a simulation in the context of a 
ACDM model with global parameters given in Ref. [9], 
with the same meatball mass spectrum which we used 
with the EdS model. Interestingly, including inhomo- 
geneities to the concordance model makes the fit worse 
and suggests, as discussed in Ref. [14], a smaller value of 
Qa- Our model remains worse than the ACDM model 
even afer inclusion of inhomogeneities, but not nearly as 
dramatically so as the homogenous EdS model. 

Of course our setup is but a simple toy model. How- 



ever, we believe that our results make it explicitly clear 
that different effects of inhomogeneities can pull into the 
same direction, making the PBS depart from the GBS. 
Moreover, we believe that the two effects of the voids 
considered here are not the end of the story; to conclude 
on the viability of inhomogeneous universe models as the 
possible explanation of the apparent acceleration, it is 
crucial to consider all possible effects in as realistic a 
model as possible. 
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